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Abstract. Microquasars, X-ray binaries displaying 
relativistic jets driven by accretion onto a compact 
object, are some of the most efficient accelerators in 
the Galaxy. Theoretical models predict Very High 
Energy (VHE) emission at the base of the jet 
where particles are accelerated to multi-TeV energies. 
This emission could be detected by present IACTs. 
The MAGIC telescope observed the microquasars 
GRS 1915+105, Cyg X-3, Cyg X-l and SS433 for 
- 150 hours in total from 2005 to 2008. We triggered 
our observations by using multi wavelength informa- 
tion through radio flaring alerts with the RATAN 
telescope as well as by ensuring the low/hard state 
of the source through RXTE/ASM and Swift/BAT 
monitoring data. We report on the upper limits on 
steady and variable emission from these sources over 
this long period. 

Keywords: gamma rays: observations — X-rays: 
binaries — stars: individual (GRS 1915, Cyg X-l, 
Cyg X-3, SS433) Radiation mechanisms: non-thermal 

I. Introduction 

Microquasars (MQ) are X-ray binary systems contain- 
ing neutron stars or black holes that display relativistic 
jets HI. Currently ~ 16 of these sources are known 
in our Galaxy Q. Since their discovery, MQs have 
attracted the attention of the high energy astrophysics 
community because they are scaled-down versions of 
quasars, and allow the study of jet emission and forma- 
tion on much shorter timescales. An evidence of TeV 
emission has been obtained in the stellar-mass black 
hole MQ Cyg X-l 0. AMQ scenario has also 
been suggested for the recently detected TeV sources 
LS 5039 and LS I +61 303 Q though they could 
otherwise contain a young non-accreting pulsar and be 
powered by its rotation [8] [9|, as is the case for PSR 
B1259 -63 OH. The detection of new MQs in the TeV 
domain is therefore an important target for the present 
VHE gamma-ray observatories. 

We present the MAGIC observations of four well- 
established galactic MQs: GRS 1915+105, Cyg X-l, 
Cyg X-3 and SS 433. These sources present unique prop- 
erties concerning the non-thermal activity, the kinetic 



power present in their jets and the surrounding radia- 
tion/matter fields where gamma-rays could be generated 
both in leptonic and hadronic scenarios. The observa- 
tions span four cycles of operations of the telescope 
from 2005 to 2008, resulting in a variety of physical 
conditions (regarding both orbital phases and spectral 
states) in which VHE emission could be produced. 

In Sect. 2, we refer the MAGIC telescope perfor- 
mance, the data reduction method and analysis technique 
applied. In Sect. 3 we present the data sample and the 
trigger strategy actually used in our observations. We 
report the resulting upper limits to the TeV emission 
for each source, both for the whole sample and in a 
phase-folded basis, in Sect. 4. Day-by-day analysis for 
transient emissions is also shown. Sect. 5 is devoted to 
the discussion of the obtained results. 

II. The MAGIC Telescope and 
Analysis Technique 

The MAGIC telescope ifTTI is a new generation Imag- 
ing Atmospheric Cherenkov Telescope at La Palma, 
Canary Islands, Spain (28.3°N, 17.8°W, 2240 m a.s.l.), 
successfully operating since the beginning of 2004. The 
parabolically-shaped reflector, with its total mirror area 
of 236 m 2 , allows MAGIC to sample a part of the 
Cherenkov light pool and focus it onto a multi-pixel 
camera, composed of 576 photomultipliers (PMTs). 

The output pulses of PMTs are converted into optical 
signals, transmitted via optical fibers and digitized by 
a Flash Analogic to Digital Converter (FADC) Sys- 
tem. In the first years of the MAGIC operations, 300 
MSample/s FADCs (Siegen-FADCs) were used and it 
was upgraded to 2 GSample/s ultra-fast FADCs (Mux- 
FADCs) in February 2007, reducing contamination of 
night sky background and improving the reconstruction 
of the timing characteristics of the recorded images. 

The data analysis was performed using the standard 
MAGIC analysis software. After the calibration lfl2l 
and the image cleaning [13], the shower images were 
parametrized by the so-called Hillas image parameters 
03). In the case of data taken with the Mux-FADCs, 
two additional parameters, namely the time gradient 
along the main shower axis and the time spread of 
the shower pixels, were also computed [13|. Hadronic 
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background suppression was achieved using the Random 
Forest (RF) method [15] where for each event the so- 
called "hadronness" was computed based on the Hillas 
and the time parameters. "Hadronness" is a measure of 
the probability that the event is not gamma like. The 
RF method was also used for the energy estimation. 
The calculation of the number of excess events was 
performed in two different ways in the data before and 
after the installation of the Mux-FADCs: in Siegen- 
FADC data, DISP method COD was used, while Alpha 
method lfl3ll was used in Mux-FADC data so that 
the timing parameters could have been used with less 
systematic uncertainties. Data taken with significantly 
lower or higher trigger rate and with zenith angle larger 
than 50° were rejected from the analysis. Images with 
less than 200 ph.e. were discarded in order to assure the 
better background rejection and to avoid a systematic 
error in estimation of the background. The sensitivity of 
the telescope is estimated to be ~ 2% and 1.6% C.U. 
for 5a detection in 50 hours of low zenith observations, 
with Siegen-FADCs and Mux-FADC, respectively IT3l 
ED. 



III. Trigger Strategy and Data Sample 



In 2005 and 2006, the MAGIC telescope observed 
GRS 1915+105 and Cyg X-3 for 14 hours each. These 
observations were triggered by the alerts on the source 
flaring state at radio frequencies sent by the RATAN-600 
telescope. On the other hand Cyg X-l was monitored 
between July and November 2006 for 46 hours and 
an evidence of signal at the level of 4.0 a (pre-trial) 
was observed on September 24th, 2006 Q, when the 
source was in a low/hard state and in coincidence with 
a hard X-ray flare. Following this promising result, in 
2007, we moved to a monitoring approach when the 
sources stayed in the low/hard state, characterized by the 
presence of steady relativistic jets where VHE emission 
could be produced. The spectral state was defined by 
using public soft and hard X-ray data from All Sky 
Monitor (ASM) iflSl . hosted by the Rossi X-ray Timing 
Explorer (RXTE) satellite, and from the Burst Alert 
Telescope (BAT) lfl9l on board of Swift, respectively. 
Table U shows the used trigger constraints on ASM [2- 
12 keV] and BAT [15-50 keV] fluxes for each of the 
observed sources. Cyg X-3 was observed in 2007 for 
28 hours; GRS 1915+105 for 14 hours in June 2007, 
and Cyg X-l for 21 hours between July and November 
2007. SS 433, was observed during the second half of 
August 2008 for 15 hours. The precession of the system 
can provoke a periodic (P* ~ 162 d) attenuation of the 
VHE gamma-rays produced near the compact object. We 
observed the source following a predicted minimum of 
the opacity [20|, avoiding also the eclipsing companion 
which cover the jet inner regions each ~ 13 d. 



TABLE I 

Trigger criteria for the microquasar observations in 
2007 





BAT daily 


BAT orbital 


ASM 




[counts/s] 


[counts/s] 


[counts/s] 


Cyg X-l 


/ 


> 0.2 


< 200 


Cyg X-3 


> 0.05 


/ 


< 200 


GRS 1915+105 


> 0.08 


/ 


< 1000 



IV. Search for steady and transient VHE 

EMISSIONS 

No signal of steady emission has been found from any 
of the four sources; in Table [TT] we report the integral 
upper limits on their VHE emission. Their values vary 
between 0.6 to 2.6% C.U, and are calculated for E 
> 250 GeV at a 95% confidence level assuming a 
power law spectrum with a spectral index of —2.6. The 
differential upper limits are shown in TableUTTl To search 
for transient emission, data have been analysed also 
night by night. There are 28, 21, 16, and 5 nights of 
observations with good conditions for Cyg X-l, Cyg 
X-3, GRS 1915+105 and SS 433, respectively. The 
distribution of the significance of the 70 observations 
is shown in Fig. [T] There is one observation with 4.0 
a excess, which is Cyg X-l observation mentioned 
above. Apart from it, the significances follow the normal 
Gaussian showing no signal evidence. The upper limit 
on VHE gamma-ray flux > 250 GeV for individual night 
is - 10% C.U. 

Phase-folded analysis has also been performed (see 
Table IIVI ). to search for possible variations of the 
gamma-ray emission at different orbital phases. No 
significant excess was found at any phase. 

V. Discussion 

VHE gamma-rays are expected to be produced in a MQ 
scenario. The emission should be naturally produced 
in the interaction between jet accelerated particles and 
the matter and radiation fields in the binary system 
and surrounding medium (see, e.g. ETIl — J24)). The 
variability and the spectral properties in these sources 
can be however very complex, depending on multiple 
conditions and parameters largely unconstrained. To 
estimate the TeV emission from these systems one 
has to take into account their intrinsic capability to 
accelerate particles to multi-TeV energies, the emission 
outcomes through different radiation mechanisms and 
the absorption of the emerging VHE flux, which can be 
strong enough to attenuate most of the putative signal at 
a given orbital phase. 

The observations presented in this work have been 
performed following a trigger strategy intended to ensure 
the low/hard state of the sources, where jets are thought 
to be present ll25l ll26ll . Furthermore, our results span 
a wide range of the phases in which each source can 
be found during its orbital motion. The phase-folded 
reported upper limits can then be used to constrain the 
VHE emission in conditions where absorption should be 



PROCEEDINGS OF THE 31 st ICRC, LODZ 2009 



3 



TABLE II 

Steady emission flux integral upper limits above 250 
GeV 





Eff. 






U.L. 




Time 






[10- 12 cm- 2 s^ 1 ] 




M 






/ [C.U.] 


Cyg X-l 


36.2 


-49.6 ±51.7 


-0.96 


1.24 / 0.75% 


Cyg X-3 


32.8 


-56.4 ± 47.1 


-1.19 


1.01 / 0.61% 


GRS 1915 


22.3 


-94.0 ± 52.3 


-1.87 


1.17 / 0.71% 


SS 433 


6.6 


-9.7 ± 27.9 


-0.35 


4.30 / 2.61% 



Significance Distribution 




Fig. 1. Significance distribution of 70 night observations of Cyg X-l, 
Cyg X-3, GRS 1915+105 and SS 433. A Gaussian distribution with 
the mean of and the standard deviation of 1 for 70 samples is also 
shown. The Cyg X-l TeV flare evidence appears at a 4.0 a level. 

the minimum allowed in each case, though the precise 
attenuation factor will depend on the intensity of the 
companion and disk photon fields, the precise location 
of the emitter with respect to them and the interaction 
angle between the gamma-ray and the target photons 
(see, for instance, [27 1 [28]). 

In a leptonic framework, Inverse Compton (IC) pro- 
cess might produce detectable VHE fluxes when scat- 
tering off external photon fields permeating the jet and 
surroundings. In addition, one should consider the con- 
tributions from IC on self-generated synchrotron photon 
fields as well as relativistic Bremsstrahlung, which could 
significantly enhance the total VHE output. 

In a hadronic scenario, gamma-ray fluxes could be 
generated through the decay of 7r° produced in the 
interactions of relativistic protons in the jet against 
thermal ions from the companion wind or within the jet 
itself. Fluxes up to few % of the jet kinetic luminosity 
could be radiated at VHE (H ED (29). 

The integral flux upper limits 0e >250 for the central 
regions of the observed sources reported in Tables [TT] 
and |IV] could point to some physical mechanism(s) by 
which the putative VHE fluxes would be suppressed. 
Synchrotron losses in the magnetic field present in the 
jets, for instance, can play a major role in the final 
spectral energy distribution of MQs. A high enough 
magnetic equipartition fraction ?7 cq = it ma g/ukin, where 



"mag and Ukin are the magnetic and kinetic energy 
densities in the jet, respectively, would make this channel 
dominate the radiative losses and the gamma-ray fluxes 
produced through IC processes would be reduced. Con- 
sidering only the companion's radiation energy density 
w* = L^/AttcPc , where d? = (a 2 rb + Zq), being L* 
the companion luminosity, a or i, the orbital separation 
at a certain phase and zq the emitter location along 
the jet axis, the non-detection of steady gamma-ray 
emission could imply u mag > u t , leading to values of 
j hq > 1.9 x lCr 3 x (L,39/iki„3 7 )(V'/0.5 o ) 2 [z 2 /(z 2 + 
)], where ip is the jet semiaperture angle, L+^q — 



7c<| 
J 2 

'orb 



"> 39 erg s \ and L 



kin 37 



= L 



kin 



/10 37 erg s" 1 is 



Wio d 

the jet kinetic power. 

In addition, the gamma-ray luminosities depend lin- 
early on q acC ch the fraction of Lkm transferred to an 
accelerated e/p plasma, which is poorly constrained 
under a theoretical point of view. It could be much 
less efficient than expected, though it seems unlikely 
since acceleration is actually required to explain the 
persistent leptonic non-thermal emission at lower radio 
energies and probably at hard X-ray energies. Interest- 
ingly, Reynoso et al. [20| have recently estimated the 
gamma-ray production in the hadronic jets of SS 433. 
By using the HEGRA upper limit on the VHE flux from 
this system [30], q acC c\ for protons is constrained to be 
< 3 x 10~ 4 . These values assume however that HEGRA 
observations took place during unknown precessional 
phases VP. In order to account for the absorption of 
the gamma-ray photons, a ^-averaged flux was used to 
derive the upper limit on q a ccoi- The current MAGIC 
observations were performed during the lowest absorp- 
tion phases. By comparing the proton-proton expected 
gamma-ray fluxes with our upper limits, a more stringent 
constraint of the hadronic q acco i < 7.4 x 10~ 5 can be 
derived. 

In addition to the steady emission, short-lived ejec- 
tions can also occur in the transitions between the 
low/hard to the high/soft MQ states, although the flaring 
behaviour is not limited to the state changes in sources 
like GRS 1915 +105 or Cyg X-3, and is quite unpre- 
dictable. The kinetic power expelled in such episodes 
can be up to several times 10 39 erg s _1 [31], producing 
intense radiation fluxes at all wavelengths. The appear- 
ance of outbursts at radio/IR wavelengths can be related 
to an increase of the activity also at VHE, although 
the time-delay between them could range from hours 
to days Oil . The hard X-ray /TeV correlation should 
provide more tightening constraints, since the emission 
could be produced by the same particle population or 
with similar timescale. However, the multifrequency 
correlation could be present only at some stages of 
the flare 0, rendering the triggering of TeV outburst 
detections a quite difficult task. The non-detection of any 
transient signal could imply the high absorption present 
in the inner regions of these systems. The gamma- 
ray attenuation with the stellar and accretion disk pho- 
ton fields could nonetheless distribute the luminosity 



4 



T. Y. SAITO et al. MICROQUASAR OBSERVATIONS WITH MAGIC 



TABLE III 

Differential upper limits for the steady VHE emission. 



Energy [GeV] 


100 - 178 


178 - 316 


316 - 562 


562 - 1000 


1000 - 1780 


1780 - 3160 


U.L. 

[cm- 2 TeV- 1 s- 1 ] 


Cyg X-l 


1.73 x 10~ lu 


2.75 x lO" 11 


3.95 x 


1.97 x 10~ 12 


4.22 X 10" i;1 


8.97 x 10~ 14 


Cyg X-3 


6.54 x lO" 1 " 


4.64 x 10" 11 


5.15 x 10 - " 


1.06 x lO" 1 ^ 


4.58 x lO" 1 '' 


1.07 x 10- 1J 


GRS 1915 


2.49 x lO" 1 " 


3.03 x lO" 11 


4.57 x 10 -1 * 


1.84 x lO" 1 ^ 


4.57 x 10" 1,5 


1.14 x 10~ 1J 


SS 433 


6.30 x lO" 1 " 


6.70 x 10" 11 


7.60 x 10 - " 


3.50 x 10"" 


1.90 x 10~^ 


3.40 x 10- 1J 



TABLE IV 

Phase-folded integral upper limits above 250 GeV 



Phase 1 


0.0-0.1 


0.1-0.2 


0.2-0.3 


0.3-0.4 


0.4-0.5 


0.5-0.6 


0.6-0.7 


0.7-0.8 


0.8-0.9 


0.9-1.0 




Effective Time [h] 


4.9 


1.8 


3.2 


0.5 


3.9 


6.6 


2.6 


2.8 


4.0 


5.9 


Cyg X-l 




0.4 


-0.6 


-0.3 


0.5 


-1.0 


0.0 


-2.9 


-0.6 


-1.3 


1.4 


Period: 5.6 d 


U.L. [lO-^cm^s-*] 


7.3 


8.0 


7.3 


30 


3.7 


5.6 


2.6 


5.6 


3.0 


10 




Effective Time [h] 


3.1 


3.0 


2.8 


3.1 


3.0 


3.6 


4.5 


3.5 


3.5 


2.5 


Cyg X-3 


a 


0.1 


-1.0 


-0.2 


-1.2 


-1.4 


-0.7 


1.3 


1.1 


0.6 


-0.4 


Period: 0.2 d 


U.L. [lO-^cm-^s- 1 ] 


7.1 


3.7 


6.0 


3.6 


3.4 


5.0 


13 


13 


11 


5.8 




Effective Time [h] 


2.0 


6.0 


0.0 


0.0 


0.0 


5.1 


4.6 


3.9 


0.1 


0.5 


GRS 1915 


a 


-0.5 


-1.0 








-1.9 


-0.2 


-0.9 


-0.9 


1.0 


Period: 33.5 d 


U.L. [lO-^cm^s- 1 ] 


10 


4.7 








2.1 


4.8 


3.6 


41 


43 




Effective Time [h] 


0.98 


0.0 


0.0 


0.0 


0.0 


1.4 


0.0 


0.45 


2.2 


1.5 


SS 433 


a 


-1.5 










1.2 




-0.24 


0.1 


-0.68 


Period: 13.1 d 


U.L. [lO-^ran-^s -1 ] 


6.5 










21 




18 


9.8 


7.7 



to lower gamma-ray energies through electromagnetic 
cascades ll23l ll27l . This would increase the fluxes at 
ranges where observatories like Fermi operate, providing 
the most simultaneous triggering of the putative VHE 
gamma-ray production in these systems. 

To conclude, steady emission from GRS 1915 +105, 
Cyg X-l, Cyg X-3 and SS 433 remains undetected 
at VHE, although an evidence of TeV flare has been 
found in Cyg X-l. A complete, global description of 
the gamma-ray emission expected accounting both for 
leptonic and hadronic processes in the different sources 
and absorption conditions at each orbital phase in which 
they have been observed is beyond the scope of this 
work. Nevertheless, some relevant theoretical predictions 
can already be tested, while further observations at 
VHE energies with present ground-based Cherenkov 
telescopes and at lower gamma-ray energies with space 
missions like Fermi, will provide new information on 
the physics present in these galactic VHE laboratories. 
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